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Simple generic aspects of nuclear pairing in homogeneous medium as well as in finite nuclei are 
discussed, ft is argued that low-energy nuclear structure is not sensitive enough to resolve fine details 
of nuclear nucleon-nucleon (NN) interaction in general and pairing NN interaction in particular 
what allows for regularization of the ultraviolet (high-momentum) divergences and a consistent 
formulation of effective superfluid local theory. Some aspects of (dis)entanglement of pairing with 
various other effects as well as forefront ideas concerning isoscalar pairing are also briefly discussed. 
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I. INTRODUCTION 



The pioneering works of Bohr, Mottelson and Pines Belyayev and, in particular, the uhimate success of 
early "large-scale" BCS calculations by Nilsson and Prior Q explaining simultaneously odd-even mass staggering 
' and moments of inertia neatly settled pairing just at the heart of nuclear physics. The success of Nilsson-plus-BCS 
approach using a simple constant matrix element pairing interaction confirms the absolute dominance of the s-wave 
^ (monopole) part of the in-medium effective particle-particle (p-p) NN interaction at the Fermi energy. There is, 
\^ ' however, yet another very general argument pointing toward the simplicity of nuclear pairing. Indeed, since pairing 

• modifies the nucleonic motion essentially only in the closest vicinity of the Fermi energy, 

it gives rise to uncertainty in momentum space, dp A/vp, what translates to uncertainty in coordinate space of the 

order of ^ ~ 1mc^}A ^ ^"^^ "^^^ quantity f , which is known as the coherence length, defines the spatial extension of 
' the nucleonic Cooper pair. Note, that [in fact, due to Ep ^ A] the value of ^ exceeds by far the typical interaction 
[ range S, ^ ^ -j^ what is known as weak coupling limit. 
<^ ' The above argumentation holds also for finite nuclei where ^ ~ 2R 3> ro [R denotes nuclear radius] i.e. nucleonic 
[ Cooper pairs are spatially very extended objects. They are therefore not bosons, they are overlapping pairs which, 
'"^ I for example, cannot form Bose-Einstein condensate. Since nuclear pairing is characterized by a small parameter 
^ ' e = ^ ^ 1 it should posses an intrinsic simplicity. In particular, it should be fairly insensitive to fine details of the 
' NN interaction or, in alternative words, should be well described by a local theory. In fact, since Tq <SC 2R, arguments 
^ ' speaking in favor of a local approximation can be extended over to the particle-hole (p-h) channel as well, see Q and 
refs. quoted therein. 

The paper is organized as follows. In Sect.|n]l shall discuss the independence of ^Sq pairing with respect to details 
^ , of the NN interaction in homogeneous medium. In Sect. lIIIl I shall present possible regularization schemes leading to a 
- - I cutoff parameter independent superfluid local density approximation. In Sect. |lVI I shall brieflv discuss the perplexing 
problem of (dis)entanglement between pairing and various other effects. Finally, in Sect. 01 shall brieffy overview 
current ideas regarding isoscalar pairing. 



II. THE ^So PAIRING GAP IN INFINITE HOMOGENEOUS MEDIUM 

Several bare NN potentials are available nowadays that fit the two-body scattering data with very high precision. 
In spite of substantial differences among them they all predict essentially the same values for the ^Sq pair-gap over 
rather wide range of Fermi momenta kp, see examples in Fig. ^ It appears that the detailed NN interaction is not 
needed at all to determine the ^Sq gap, as nicely demonstrated in Ref. I^. The decisive point is that the ^Sq NN 
scattering is characterized by a large negative scattering length indicating the presence of a nearly bound resonant 
state at zero scattering energy. Around this low-energy pole the NN interaction can be well approximated by a 
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FIG. 1: ^So pair-gap caculated using CD-Bonn and Nijmegen bare NN interactions and free spectrum. Dotted line shows 
calculations using phase-shift approximation. Dash-dotted line indicates results obtained by using finite range approximation 
for phase shifts. Taken from Ref. 0]. 



separable interaction v{k, k') ~ Xv{k)v{k') for which the so called inverse scattering problem can be solved. It means 
that both v{k) as well as A(fci?) are in fact fully determined by means of the phase shifts 5{k). Although, at first 
glance, this approximation seems to be valid (?) only at low energies and requires {ii) the knowledge of 5{k) at, in 
principle, all energies it appears to work surprisingly well up to kp ^ 1.4 fm~^. Note however, that the effective range 
approximation to the phase shifts works well only at low-densities up to kp ~ 0.6 fm~^, see Fig.^ 

The inclusion of in-medium polarization and screening corrections appear to be extremely difficult. So far no 
consensus has been reached how to consistently compute these corrections, see Refs. ii and refs. therein. However, 
instead of deriving in- medium gap equation lS.{kp) from the bare NN interaction one can attack the problem starting 
directly from an effective interaction like Gogny or local density dependent delta interaction (DDDI). Indeed, as 
shown by Garrido et al. ^ the pair gap calculated using the Gogny force with parameter set DIS fitted directly to 
finite nuclei follows rather closely the pair-gap calculated using the bare Paris-force with Brueckner-Hartree-Fock 
spectrum. The difference between the gaps increases with increasing density reflecting, most likely, an enhancement 
due to the averaging over low-density skin region which is implicitly taken into account through the parameters of the 
effective Gogny interaction. In Ref. |3] it was demonstrated that DDDI is also flexible enough to follow the Paris-force 
(and Gogny) results, however, only after careful adjustment of the cutoff parameter. A need for cutoff within the local 
pairing theory is neither satisfactory nor unambiguous and spoils to large extent its simplicity. It is interesting to 
observe, however, that local interactions can be safely used in p-h channel [Skyrme force] without any artificial cutoff 
and that in fact in this channel it is equivalent, at least according to effective theory principles, to the finite-range 
Gogny interaction. 



III. TOWARD LOCAL SUPERFLUID EFFECTIVE THEORY. 



A guiding principle underlying any effective theory aiming to describe the low-energy limit of a deeper, more 
fundamental theory can be formulated in the following way: The low-energy (infrared) phenomena are not sensitive 
enough to resolve high-energy (ultraviolet) dynamics. It means that short-range (high-momentum) dynamics can 
be removed from the theory and replaced (renormalized) by a few local corrections. In momentum space one can 
therefore expand a short-range (SR) interaction as: 

vsR{q^) -9 + 92q^ + 94q* ■■■ , (1) 

expressing it formally by means of few constants 9,92,94, ■ ■ ■ which need to be carefully readjusted to a selected set 
of low-energy data. In coordinate representation interaction Eq. ^ can be rewritten as: 

vsR{r) « ca'^Sa{r) 

+ dia^V^Sair) + d2a^V5a{r)V 

+ ... 

+ h,a-+'V'^Sa{r) + ... , (2) 
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where Sa{r) denotes an arbitrary model of the Dirac delta function while c, di, ^2, ^i, • ■ ■ denote empirically adjustable 
constants. In nuclear structure, due to the non-singular and very short range nature of the effective interaction, one 
can in fact use the strict limit \ima-,Q Sa{r) = 5{r) which corresponds to the well-known Skyrme interaction [first 
three terms in Eq. Q]. On the other hand by retaining in Eq. Q only the first term modeled by a sum of attractive 
and repulsive Gaussians of different ranges one obtains the well-known finite-range Gogny force. Of course both forces 
must be augmented by the density dependent and spin-orbit terms as well as by space-, spin-, and isospin-exchange 
terms. The Skyrme and Gogny forces are therefore two realizations [among infinitely many] of effective nuclear 
interactions which are used in practical nuclear structure calculations. 

The Gogny force can be unambiguously used also in the p-p channel. Indeed, the finite-range (tq ^ 1 fm) auto- 
matically discriminates states above Ec ^ pl/2mr ~ /mro ~ 40MeV [m^. = m/2 is reduced mass] since VoPc ~ ^■ 
The Gogny interaction was proved to be very successful in numerous practical applications. In particular, the average 
Hartree-Fock-Bogolyubov (HFB) pair-gaps calculated using DIS Gogny follows very accurately empirical three-point 
OES A(iV) = {-l)^[B{N - 1) + B{N -f 1) - 2B{N)]/2 calculated for odd-A^, see Ref. [13 and Fig.H 
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FIG. 2: Mean neutron pair gaps caclulated using HFB DIS Gogny method (open dots) in Ref. [13] in comparison to empirical 
three-point OES data A(orfd - iV) (filled dots). 



Local delta-type or DDDI type pairing interactions are successfully used particularly in conjunction with Skyrme 
interaction. These applications explicitly invoke a cutoff parameter to avoid divergences. It appears, however, that 
the divergence can be rather easily identified and subsequently regularized and a cutoff free local superfluid theory 
can be constructed. 

To illustrate a possible regularization schemes let us consider first the effective theory approach of Papenbrock and 
Bertsch ^3]. In the contact-force approximation which takes into account only first term in Eq. the BCS gap 
equation takes the following form: 



2(27r) 



(Pk 



(3) 



where V and A denote volume and chemical potential, respectively. The integral in Eq. is ultraviolet (high- 
momentum) divergent. It appears, however, that the equation for scattering length, a, is also divergent within 
contact approximation: 

(4) 
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and that the divergences are of the same type. Hence relation (0J can be used as a counter-term to regularize the gap 
equation: 



m 
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^{Sk - A)2 + A2 



(5) 



This is a very elegant example of regularization connecting the pair-gap (and contact-force strength g) directly 
to the free two-particle scattering length. However, the formalism applies only to dilute homogeneous medium. 
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FIG. 3: Experimental OES (dots, solid lines) versus calculated (triangles) OES (upper part) and calculated value of mean 
proton LN gap Aljv — A2 (lower part). Filled (open) symbols denote GS (HSI) values, respectively. Gray triangles include 
theoretical GS values of OES corrected by residual pn interaction. 



Finite range corrections [i.e. higher order expansion terms in Eq. are difficult to handle but can be, at least in 
principle, systematically implemented and regularized order by order. Judging from Fig. ^ the lowest order finite 
range corrections are expected to extend the validity of this scheme till kp ^ 0.6 fm"^. 

The problem of ultraviolet type divergence in anomalous density matrix persist also in finite nuclei: 

J^(r-i,r2) = V<(ri)u,(r2) - I (6) 

Here the situation seems to be even more complex since (i) realistic single-particle (sp) spectra must be used right 
from the beginning and (ii) it is not at all obvious what physical quantities need to be used in order to regularize 
divergent terms. The appropriate regularization scheme was proposed recently by Bulgac and Yu |l2l Il3j. Their 
scheme is built upon the local density approximation (LDA) i.e. takes automatically into account the dominant 

p-h channel. The idea is to introduce cutoff {Ec = -^^^7-) dependent counter-terms leading to standard local HFB 
formalism with cutoff parameters but with a gap equation dependent on the effective running coupling constant: 

^c{r) = v*{r)ui{r) , 

Ei>0 

Mr) = -geff{r)u,{r), 
1 1 ■m{r)kc{r) ( 

Introducing a running coupling constant implies that the cutoff dependence is only formal and disappears for suffi- 
ciently large Ec 0,E3- This cutoff free superfluid LDA (SLDA) approach is now in phase of extensive tests [T^IT^ . 



kpir) 



2fcc(r) fc,(r) ~ kpir) 



(7) 

(8) 
(9) 



IV. ENTANGLEMENT 



Pairing gaps are not directly accessible in experiment. Hence, various indirect methods must be applied to extract 
information about them. The major difficulty is that all these indirect methods entangle pairing with various effects 
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FIG. 4: Calculated energy differences AEg^p — AEth between terminating states, AE = E{d'^^^f"^^) — -B(/"/2)' i^i A~50 mass 
region relative to experimental data. Open symbols denote state of the art SM results (shifted by 475keV) while filled symbols 
denote the SHF calculations using the SkO parameterization with spin-orbit term reduced by 5%. Gray dots denote the SHE 
results in N=Z nuclei corrected phenomenologically for isospin breaking effect, see insert and Ref. [2^ for more details. 



including shape and shape-polarization effects, sp splitting, time-odd fields, or beyond mean-field residual interaction 
effects making life rather perplexing and, in fact, introducing in a natural way uncertainties into our knowledge of 
nuclear pairing. It is therefore desirable to hunt for simple physical situations or phenomena where at least some of 
these contaminants are either decoupled or can be relatively well controlled. Superdeformation (SD) is one of the 
most prominent examples of such a phenomenon. Indeed, stability of nuclear shape along the SD band allows to study 
pairing correlations from the static to dynamic regime. Let us recall that such concepts and techniques like double- 
stretched quadrupole pairing "Ts'l, the surface-active DDDI "T^, the Lipkin-No garni (LN) number-projection Tt] were 
applied for the first time in a systematic way in SD bands in Hg-Pb nuclei |l5l llSl Il9l |20j . Afterwords these methods 
became standard in large-scale calculations in high-spin physics. 

High-spin isomers (HSI) open yet another and so far unexplored venue to study pair correlations and blocking 
phenomena [21. 22j . Thanks to their structural simplicity both configuration and shape can be kept rather well under 
control. In turn, shape and pairing polarization due to blocking can be studied in detail and traditional average gap 
method '23^ used to determine pairing strength, Gmn, can be re-examined. In particular, it was found that inclusion 
of these polarization effects requires ~10% larger G as compared to Gmn |^ to reproduce experimental OES. 

Recently the HSI have been systematically observed in A''=83 nuclei with 60<Z<67 This unique data set 

enables to study for the first time OES both at the ground states (GS) as well as at high-spins. The most striking 
feature of this data set is the almost constant excitation energy of the HSI what implies that /S.gs{Z) « Ahsi{Z), see 
Fig. 121 Conventional interpretation of this result in terms of pairing-gap suggests the lack of blocking phenomenon. 
The pair-gaps A^^v — A2 calculated using diabatic Strutinsky type method involving self-consistent blocking and LN 
particle- number projection |22l |24| shown in lower part of Fig. 13 clearly show that this conventional interpretation is 
oversimplified. Contributions to OES from pairing (blocking), sp-proton energy splitting and residual proton-neutron 
(pn) interaction must all be taken into account to reproduce experimental data in a satisfactory way, see upper part 
of Fig. |51 Even then, contribution due to blocking is clearly too strong and requires revisiting. 



V. HUNTING FOR FINGERPRINTS OF ISOSCALAR PN PAIRING COLLECTIVITY. 

The existence of isoscalar (t = 0) pairing phase in atomic nuclei is still a fascinating open issue. Various phenomena 
are discussed in this context including an onset of < = pairing driven by nuclear rotation. However, neither 
the suggested shifts in ground-band S-band crossing frequency nor substantial changes in moment of inertia were 
convincingly confirmed by experiment. On the contrary standard calculations seem to work reasonably good in 
particular in the most promising area of heavy deformed A~80, N ^ Z nuclei 113, H^- There are, however, new data 
indicating rather unusual band crossing phenomena in e.g. '''^Kr |29l |. Indeed, standard mean-field calculations p9l | can 
reproduce the ''■^Kr data assuming two different structures below [one-quasi-particle (QP) band built upon negative 
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parity Nilsson-level originating from v{f^j2-P'3/2) subshell] and above back-bending where odd-neutron occupy positive- 
parity Nilsson state originating from 1^59/2 subshell while protons form 2QP structure involving 7r(/5/2-P3/2) ® 1^99/2- 
Strong E2 transitions connecting these structures (i) cannot be explained within standard mean- field calculations 
and {ii) indicate unusually strong configuration mixing. Whether or not this configuration mixing can be accounted 
for within mean-field approximation invoking t = pn-pairing i.e. novel spontaneous symmetry breaking mechanism 
resulting in pn- mixing is under study [3(1 |. 

Also calculated energy differences Ai? = -£'('^3/ 2/7/2^) ~^(/7/2) between band-terminating states in A~50 seem to 
provide new evidence for t = pairing [2^. As shown in Fig. ^ the values of AE^xp — AEth calculated using state of 
the art shell-model (SM) and Skyrme-Hartree-Fock (SHF) method follow, up to a constant offset of ^475 keV, very 
closely each other. The quantitative difference between the accuracy of theoretical predictions 'm N ^ Z as compared 
to N = Z nuclei seen both in the SM as well as in the SHF calculations suggests, most likely, an enhanced t — pair 
scattering from sd to fp shell which is beyond present shell-model space. Similar effect was suggested in Ref. |3ll |. 

One of the most promising signal of t = pairing comes from binding energies 'm N = Z nuclei, where the problem 
of the Wigner energy (WE) [extra binding energy] is known to plague mean-field masses, see Ref. js^ and refs. 
therein. It is relatively well established from nuclear shell-model studies that the WE is predominantly due to i = 
interaction [ssLls^ l. The SM indicates, however, that the structure of the WE is very complex and is not dominated 
neither byL = 0,S' — 1,T = [s^ nor by J = 1,T = [sj isoscalar pn-pairs. Within mean-field model, which 
is entirely based upon the concept of spontaneous symmetry breaking, the definition oi t ~ Q pairing in terms of 
symmetry conserving pairing interaction consisting only L = Q, S — l,r = (or J — \,T = Q) pairs is not at all 
justified. Although the form of effective i = pn-pairing interaction appropriate for mean-field calculations is an open 
issue it is rather well established qualitatively that mean-field model augmented by t = pairing is capable to heal the 
mass-defect problem around N ^ Z line ,^6, 37, 38J. Quantitative estimate oit — Q pairing within mean-field requires, 
however, a reliable evaluation of corrections due to isospin-symmetry restoration which go beyond mean-field. Within 
random-phase-approximation these corrections roughly restore linear term in the nuclear symmetry energy (NSE) 
giving rise to ^ T{T + 1) dependence of the NSE 39]. It appears however that at least part of the linear term (~ T) 
is incorporated already at the level of (self-consistent) mean-field. Hence, further progres s in the field is impossible 
without thorough understanding of the NSE being currently under intense studies |4Cl |4 ll |42L |43| . 
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